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Recent progress in unquenched lattice QCD simulations is reviewed with emphasis on understanding of chiral 
behavior for light quark masses. 



1. Introduction 

In the long standing challenge of the calcula- 
tion of hadron spectrum from the first-principles 
lattice QCD simulation, the inclusion of dynami- 
cal quarks, or unquenching, has been the crucial 
but computationally demanding step. Thanks to 
the rapid increase of the computer power and re- 
cent developments in the lattice formulations and 
numerical algorithms, such simulations have be- 
come feasible in the past several years. 

One of the stumbling block toward full QCD 
simulations has been the treatment of the third 
dynamical quark, i.e., the strange quark. Since 
the standard HMC algorithm pP assumes that the 
number of flavor is even, previous simulations in- 
cluding the odd number of flavors had to use a 
non-exact algorithm, such as the R algorithm [2], 
which contains a systematic error of order St'^, 
with St the molecular dynamics step size. Re- 
cently an exact HMC algorithm with polynomial 
approximation has been developed 3 and this 
difficulty has been essentially removed. 

A major problem in the dynamical fermion sim- 
ulations is that the up and down quark masses are 
much smaller than the masses we can simulate on 
the lattice. Therefore, we usually extrapolate the 
data obtained at an order of magnitude larger sea 
quark masses to the physical point assuming some 
fitting ansatz. For the pseudo-sealer meson chan- 
nel, the fitting ansatz can be justified by Chiral 
Perturbation Theory (ChPT), which provides a 
systematic expansion in terms of pion mass and 
momentum squared. However, the problem is 
that the region of applicability of ChPT is not 
known a priori, and the lattice data, for instance 
obtained by the JLQCD collaboration, do not in- 



dicate the non-analytic behavior expected from 
the next-to-leading order ChPT within simulated 
quark masses, e.g., above ms/2 for JLQCD 0]. 

Since the last year's Symposium, a number of 
unquenched lattice QCD simulations have been 
pursued to overcome the problems related to the 
chiral behavior of hadron spectrum, which I re- 
view in this talk. In Section |21 1 summarize the 
major unquenched simulations available to date 
for both 2- and 2-|-l-flavor cases. In Sectional 
discuss the algorithmic issues and physics results 
from recent 2+1-flavor simulations. Discussion on 
the chiral extrapolation, which is the main focus 
of this review, is given in Section^ 

The theoretical background of the chiral per- 
turbation theory is discussed by Bar |5l6| at this 
conference. For the unquenched calculations of 
quark masses, weak matrix elements, and heavy 
quark physics, see other reviews in |7I8I9| . 

2. Recent Unquenched QCD simulations 

Table ^ summarizes the recent large-scale un- 
quenched simulations with Nf dynamical flavors. 
The lattice spacing a, physical spatial extent of 
the lattice aL, and pseudo-scalar and vector me- 
son masses Mps and My at the unitary point 
{i.e., valence quark mass equals sea quark mass) 
are listed. 

Throughout this paper we use the following ab- 
breviations to denote the lattice actions. The 
gauge and quark action combinations are denoted 
as, PW: plaquette gauge and unimproved Wilson 
quark actions, PC(X): plaquette gauge and clover 
quark actions with the improvement coefficient 
csw determined by a method X, RC(X): RG- 
improved gauge and clover quark actions with the 
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Table 1 

Recent spectrum runs in dynamical QCD simulations. an update of |10lllj . ^ an update of \V2\ ) 



CoUab. 



Nf 



Action 



fl [fm] aL [fm] AU/Mp A'U [GcV] 



Rcf. 



MILC 

CP-PACS 
& JLQCD 



2+1 



2+1 



SZAT 



RC(N) 



0.125, 0.09 
0.1 



„ . o n 0.58-0.30 
0.49, 0.38 

1.6, 2.0 0.77-0.64 



0.60-0.25 



1.00-0.64 



CP-PACS 
CP-PACS 
UKQCD 
UKQCD 
JLQCD 
QCDSF 
& UKQCD 
RBC 
qq+q 
SPQcdR 
GRAL 
SESAM 
SESAM 
& TxL 



RC(T) 
RC(T) 
PC(N) 
PC(N) 
PC(N) 

PC(N) 

DBW2DW 
PW 
PW 
PW 
PW 

PW 



0.22-0.09 
0.22 
0.1 
0.1 
0.1 

0.1 

0.12 

0.20 
0.066 
0.13, 0.08 
0.086 



2.6-2.1 0.8-0.5 

2.6 0.61-0.35 

1.6 0.84-0.58 
1.6 0.44 

1.1- 1.8 0.80-0.60 

1.2- 2.4 



1.9 

3.2 
1.1,1.6 
1.1-2.1 

1.4 



0.092, 0.076 1.3-1.8 



0.6-0.8 

0.83-0.69 
0.83-0.57 



1.2-0.5 
0.58-0.27 

1.1- 0.6 
0.42 

1.37-0.60 

1.2- 0.64 

0.49,0.61,0.70 
0.66-0.37 

0.42-0.64 
1.00-0.64 

0.90-0.49 



|1()IHj 

CHI' 
US 

m 
m 

1221 

123 
EH 
123 



improvement coefficient determined by a method 
X, SZAT : Symanzik improved gauge and AsqTad 
KS quark actions 0(a^ )-improved at tree level, 
DBW2DW : DBW2 improved gauge and domain 
wall quark actions. The method X for clover 
quark action is either N: non-perturbatively de- 
termined, or T: tadpole estimated. 

Let us touch upon the representative full QCD 
simulations. For the Wilson-type quark action, 
the continuum extrapolation was previously at- 
tempted by the CP-PACS collaboration jlOlllj 
using data at three lattice spacings with Nf = 2 
RC(T) action, and a better agreement of hadron 
masses than in the quenched case was observed 
in the continuum limit. Ref. |18j is an update 
of this work, which extends the simulation at 
the coarsest lattice spacing (a ~ 0.22 fm) to- 
ward smaller sea quark masses corresponding to 
Mps/Mv =0.60-0.35. 

UKQCD 2^ also recently extended their previ- 
ous simulation ^2) toward smaller quark masses. 
The SPQcdR collaboration reported their prelim- 
inary results for the light hadron spectrum, light 
quark masses and renormalization constants |25| . 
The GRAL collaboration 26 has started a study 



of the finite volume effect at small quark masses 
on small to medium-sized lattices. All these are 
Nf ^ 2 simulations. 

The CP-PACS and JLQCD collaborations have 
jointly started Nf = 2 + 1 flavor simulations with 
the RC(N) action, and their first results for a 
16^ X 32 lattice appeared at the last lattice confer- 
ence jl6j . They have extend the physical volume 
to 20^ X 40 this year. 

The Nf = 2 and 2 + 1 unquenched simulations 
with the KS-type fermions can be performed at 
smaller quark masses compared to the Wilson- 
type quark actions, because the required compu- 
tational cost is much reduced. The 2+1-flavor 
simulations with the SZAT action have been car- 
ried out for several years by the MILC collabora- 
tion |13ll4ll5j . This year they reported a detailed 
analysis of the chiral and continuum extrapola- 
tions by adding data at smaller sea quark masses 
and finer lattice spacing. 

With the configuration generated by the MILC 
collaboration, the "gold-plated" hadronic observ- 
ables have been calculated by the HPQCD-MILC- 
UKQCD-Fermilab collaborations and good agree- 
ment with experimental results is observed |30| . 
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The RBC collaboration reported the calcu- 
lation of pseudo-scalar meson masses and de- 
cay constants using the unquenched domain wall 
fermion j31l22| for Nj = 2 case. 

3. Nf — 2 + 1 simulations 

3.1. Algorithmic issues 

The simulation algorithm widely used for un- 
quenched simulations is the Hybrid Monte Carlo 
algorithm 1 . Quarks are treated by the pseudo- 
fermion method .32^, with which we can treat 
even number of them dynamically. If the single- 
flavor lattice Dirac operator D has a real positive 
determinant det[-D], it is also possible to treat 
it as a probability distribution function for link 
variables. The naive application of the pseudo- 
fermion method to det[Z3] results in a complex 
effective action, however. This problem is avoid- 
able by constructing an operator S which satis- 
fies — D. The polynomial approximation and 
rational approximation can be used to construct 
such an operator j^Hl- With this method the 
HMC algorithm with the polynomial or rational 
approximated pseudo-fcrmion has been obtained. 
The algorithm is then combined with the usual 
two- flavor HMC algorithm to make an Nf = 2 + 1 
flavor HMC algorithm [3]. The efflciency of the 
algorithm is comparable to that of the HMC al- 
gorithm 34. 

For the KS-type quark actions, the fourth-root 
trick is widely used to express a single flavor of 
fermion. This amounts to taking a fourth-root of 
the determinant of the KS Dirac operator -DkSi 
which represents four flavors of fermions in the 
continuum. This trick is combined with the R 
algorithm [2] for carrying out simulations. 

An important issue with the R algorithm is 
controlling the systematic error of order St"^ at 
flnite molecular dynamics step size 5t. It is ar- 
gued that a condition 5t < m is needed [33] to 
avoid large systematic errors. A reliable estimate 
of the magnitude of error is difficult prior to the 
actual simulations. 

In the MILC simulations |13I15| 6t < 3 /2am 
is adopted with light dynamical quark mass am. 
The systematic error for some obscrvables is in- 
vestigated by comparing the results from addi- 



tional simulations with larger 6t. The additional 
runs are too short to conclude the (St)'^ behavior 
for hadronic masses. Further investigation on this 
issue should be made by increasing the statistics 
or using exact algorithms as described below. 

The polynomial approximation described for 
the Wilson-type quarks above can also be applied 

1 /4 

to Djlg . This leads to the exact algorithms, the 
polynomial HMC which uses polynomial approx- 
imant and the rational HMC which uses rational 
approximant, have been developed for the two- 
flavor case |36l35l37l38j . For more details of the 
recent development of the simulation algorithm, 
see EHl- 

In order to justify the fourth-root trick one has 
to show the existence of a local fermion kernel D 
which satisfies det[Z3Ks]^^^ = det[£']. Otherwise 
there is no guarantee that the continuum limit 
is real QCD. It is shown that a naive candidate 
D = D^/g is non-local |4(JI41| . and the existence 
of a local operator D is still an open question. 
The related issues are discussed in [42143144] . 

3.2. Nf = 2 + 1 simulation v^fith the Wilson- 
type fermions 

The CP-PACS and JLQCD collaborations have 
started a joint program to perform realistic sim- 
ulations including the dynamical strange quark. 
They performed simulations on two lattice vol- 
umes, IB'^ X 32 and 20'^ x 40, at a finite lattice spac- 
ing a ~ 0.1 fm with the RG- improved gauge ac- 
tion and non-perturbatively 0(a)-improved Wil- 
son action. In order to interpolate or extrapolate 
to the physical strange quark mass, the simula- 
tions are made at two strange quark masses corre- 
sponding to Mps,ss/^^v,ss - 0.71-0.77. For the 
light up and down quarks which are assumed de- 
generate, six (five) quark masses on the 16'' x 32 
(20"^ X 40) lattice are simulated. It covers the 
range Mps,ll/Mv,ll ^ 0.62-0.78. 

They calculated the light meson mass spec- 
trum, and up and down, and strange quark 
masses. The meson masses are calculated at the 
unitary point where the valence and sea quarks 
have equal masses. Light and strange mesons are 
constructed with Light-Light (LL), Light-Strange 
(LS) and Strange-Strange (SS) combinations of 
simulated quarks. Chiral extrapolation is made 



4 



4.0 



CM 



= 3.0 



2.0 



N^2+1 



RG + clover(NPT) 



^-input 

K-input 

AWI 

a~0.1 fm 



N^2 



RG + clover(PT) 



AWI 

a~0.1 fm 



RG + clover(PT) 



AWI 
a~0.1 fm 



160 



140 



§ 120 
II 



Is J 00 

s 



80 



60 



N^2+1 



RG + clover(NPT) . 



if -input ■ 
K-input ^ 

AWI 
a~0.1 fm 



N^2 



RG + clover(PT) 



K-input 



i>-input 

4 



AWI 
a~0.1 fm 



RG + clover(PT) 



^-input 



K-input 



AWI 
a~0.1 fm 



Figure 2. Nf dependence of light (left panel) and strange (right panel) quark masses at a ■ 
CP-PACS and JLQCD [n|. 
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Figure 1. CP-PACS- JLQCD result for the Nf 
dependence of the light meson masses at a ^ 0.1 
fm 22]. Circles represent the Nf = 2 + 1 result 
with RC(NP) [ni, while triangles show the Nf = 
2 and data with RC(TP) [III- 

by a simultaneous fit to LL, LS and SS masses 
with a polynomial fit function including up to 
quadratic terms of valence and sea quark masses. 

Figure ^ shows the deviation of meson spec- 
trum from the experiment value. Lattice data 
are shown for the 2-1-1 (filled circle), 2 (up trian- 
gle), and (down triangle) dynamical flavors at 
a '--^ 0.1 fm. The error bar shows statistical error 
only. Compared to the 2-flavor result, the 2-1-1- 
flavor simulation gives the result slightly closer 
to the experiment. Figure El shows the compari- 
son of quark masses with different Nf at a ~ 0.1 
fm. The reduction of quark mass with the inclu- 



sion of dynamical up and down quarks is signifi- 
cant. The effect of dynamical strange quark is, on 
the other hand, not conclusive at this stage, be- 
cause there may be an effect of a different choice 
of csw yielding different scaling violation among 
the data. The light quark mass could be sensitive 
to the choice of the detail of the chiral extrapola- 
tion. The effect of the chiral logarithm should be 
investigated. 

3.3. Nf = 2 + 1 simulation with the KS-type 
fermions 

The MILC collaboration has been generating 
configurations of iV/ = 2 + 1 QCD with the SZAT 
action. The HPQCD-MILC-UKQCD-Fermilab 
collaborations calculated the "gold-plated" ob- 
servables on these configurations, and found re- 
sults shown in Fig. A good agreement is seen 
between the experimental and lattice results for 
Nf = 3 for a number of quantities covering both 
the light sector and the heavy sector. 

This year the MILC collaboration published a 
detailed analysis of the continuum and chiral ex- 
trapolations 1451 15| . Since the analysis of the 
light pseudo-scalar meson channel using ChPT 
is very important, and also instructive to the 
Wilson-type quark action simulations, I discuss 
the key points of their method to obtain the chi- 
ral and continuum limit in Sec. 4. 

From this analysis, they obtain the light meson 
spectrum, the decay constants, and the low en- 
ergy constants of chiral perturbation theory |45|. 
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Figure 3. Lattice QCD results with SZAT ac- 
tion divided by experimental results for the "gold- 
plated" observables 150] . 

The light quark masses are also extracted |4(j) . 
The results for the pion and K meson decay con- 
stants are 

U = 129.5 ±0.9 ±3.5 MeV, 
fK = 156.6 ± 1.0 ±3.6 MeV, (1) 
/k/U = 1.210(4)(13), 

which are updated values of Fig.|21 The errors are 
statistical and systematic, the latter containing 
chiral and continuum extrapolations, the scale de- 
termination, and electromagnetic effects as care- 
fully investigated in [T5) . 

The strange quark mass rus, and averaged up 
and down quark mass m in the MS scheme at 
0{as) matching are given by 



MS 



MS 



m 
ms/rh 



76(0) (3) (7) (0) MeV, 

2.8(0)(1)(3)(0) MeV, 
27.4(1)(4)(0)(1), 



(2) 



where the errors are from statistics, simulations, 
perturbation theory and electromagnetic effects, 
respectively. 

The light hadron spectrum reported in |15j is 
also an update of Fig. |3| The systematic er- 
rors from chiral extrapolation are not included 
in Fig. 16 of 15.^, but estimated in Fig. 4 of [T^ 
for the nucleon mass. 

An intersting work is discussion of the candi- 
date of two-particle state for 0++(ao) and 1^ {hi) 



channel with the observation of the level crossing 
as decreasing the quark masses, although further 
clarification is still needed with smaller statistical 
error and smaller quark masses. 

The heavy hadron mass splittings are obtained 
in [30], where the NRQCD action and Fermilab 
actions are used for b and c quarks respectively 
as shown in Fig. 13 (and also Fig. 16 of )15|). 

A highlight of prediction from lattice QCD 
is the mass of meson which is not well es- 
tablished experimentally. The current experi- 
mental value is 6.4(4) GeV from CDF [32] or 
5.95(37) GeV from DO prehminaly data [3H|- Us- 
ing the MILC configurations Mb^ has been calcu- 
lated by the HPQCD-FNAL-UKQCD collabora- 
tion applying the same technique to the heavy 
hadron spectrum described above. They obtained 
Mb^ — 6.304(16) GeV. Comparions with future 
improved experiments will provide a verification 
of the technique used for the configuration gener- 
ation. 



4. Chiral extrapolation of hadronic observ- 
ables 

Consistency of lattice data with the chiral log- 
arithm has been an important issue since the lat- 
tice 2002 conference [20] . Possible reasons for not 
finding the logarithm are : (1) quark masses are 
still too heavy to apply ChPT, and/or (2) lattice 
cutoff effect distorts the chiral logarithm. Both 
effects could be important for available simula- 
tion parameters 0. Here we discuss recent re- 
sults bearing on this question. 

4.1. KS-type fermion action and SChPT 

The MILC Collaboration use the staggered chi- 
ral perturbation theory (SChPT) as a guide in 
the chiral and continuum extrapolations [451 15| . 
SChPT represents the broken SU{A x 3)l x 
S'C/(4 X 3)ij symmetry, and incorporates both 
the taste symmetry breaking effect and fourth- 
root trick. The one-loop formula for the pseudo- 
scalar masses and decay constants are obtained 
Using the data at different lat- 



tice spacing and partially quenched [H^ data sets, 
they simultaneously fit Afpg and /ps, and then 
the chiral and continuum limit can be taken at 
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the same time. They have two lattice spacings, 
a ~ 0.125 fm (coarse lattice) and ~ 0.09 fm 
(fine lattice). On the coarse lattice, five values 
of the light quark masses am' corresponding to 
the range of Mps,ll from 0.60 GeV to 0.25 GeV 
are simulated, and the strange quark mass am'^ 
is fixed at a slightly higher value than nature 
and the "X- meson" mass Mpg^LS ranges from 
0.69 GeV to 0.58 GeV. On the fine lattice, two 
values of am"s at a fixed am'^ are available to 
date, which correspond to AfpgxL — 0.45-0.32 
GeV and Afps,LS = 0.60-0.55 GeV (the prime on 
mases means that they are the mass of dynamical 
quarks used in the simulations). In the fitting the 
partially quenched points are used, which could 
help to stabilize the fit In the MILC analy- 
sis the valence light quark masses are in the range 
O.lm'^-m'^. 

The taste breaking effect appears in the mass 
splitting between Goldstone and non-Goldstone 
pions. With the AsqTad action the breaking of 
O(a^a^) is expected, which is confirmed by mea- 
suring the ratio of the mass-squared splitting of 
fine to coarse lattice (AAf^)fine/(AM^)coarsc — 
(aa^)fine/(Q!a^)coarsc = 0.372. The splitting 
also shows the 5*0(4) symmetry expected from 
SChPT 

The mass range where the (S)ChPT is valid 
is not known a priori. To see how large quark 
masses can be used in the chiral extrapolation 
they attempted the fit for different subsets of the 
data. The data subsets are: 

• subset I (94 data points) with rrix + ruy < 
0.4m(. (coarse) and rrix + my < O.bAm'^ 
(fine). 

• subset II (240 data points) with m^ +my < 
0.7m'j, (coarse) and m^+my < O.Sm'j, (fine). 

• subset III (416 data points) with mx+my < 
l.lOm'j, (coarse) and m^ + my < 1.14m'j, 
(fine), 

where m^^y) means the mass of valence quark. 
They found that the NLO formula is not sufficient 
even if for the subset I (lightest set), for which 
the heaviest meson mass is around ~ 0.5 GeV. 
At this value the NNLO contribution is expected 
to be around ~ 3.5%, which is larger than the sta- 
tistical error (0.1%-0.7%). Since the fuU NNLO 



terms for SChPT is not known, they include an- 
alytic terms in the fitting ansatz. The NNNLO 
effect is also investigated. For future analysis full 
NNLO form for SChPT is desired. 

The fit parameters include the effects of the 
scaling violation. Taste symmetry breaking terms 
could vary as 0(Q!^a^), and other parameters in- 
clude dependence on O(asa^) or 0{ala?') in the 
global fitting. Thus, they investigated the fol- 
lowing combination of fit ansatz and data sets to 
estimate various systematic errors. 

• NNLO fit on subset I. 

• NNLO fit on subset II. 

• NNNLO fit on subset III. 

where NNNLO fit is used to interpolate strange 
quark mass. 

Figure 01 shows the result of simultaneous fit- 
ting (NNNLO fit on subset III). Selected data 
points and fit curves are presented. After the 
global fitting, the parameters are extrapolated to 
their continuum values, then the ChPT formula 
at continuum limit is recovered. The infinite vol- 
ume limit is also taken. 

Existence of the chiral logarithm could also be 
examined. Both the continuum logarithmic terms 
and taste breaking logarithms are needed to get 
good fits. Although the finite volume effect is 
expected to be small with their physical volume 
[L > 2.5 fm), it is observed that the fit degrades 
when the finite volume correction is switched off. 
The good fit is highly non-trivial and rely on the 
good statistics and partial quenching. It is not 
merely a consequence of the large number of fit 
parameters (40 parameters for NNLO fit). 

The physical point is then derived by adjusting 
m! and with the ChPT formula at continuum 
limit to the corresponding mess masses and decay 
constants. They considered the electromagnetic 
and isospin violating effect when adjusting the 
physical point since their statistical precision are 
comparable to these effects. Thus the low energy 
constants, decay constants, quark masses are ob- 
tained |15I46| . 

They observed the chiral logarithm and ob- 
tained the best results in the chiral/continuum 
limit. For further improvement, study on scal- 
ing violation using three or more lattice spacing 
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is desired, with which the assumptions on the 
scaling of the coefficients in the ChPT formula 
can be verified. The partially quenched NNLO 
formula is now available in the continuum the- 
ory which would give hints to the NNLO 
SChPT calculations. 

They observed that the ChPT formula gives 
a natural expansion with respect to the quark 
mass and taste symmetry breaking terms (0(1) 
coefficients in the expansion). By a naive order 
counting NNLO contribution is ~ 10% even at 
Mps ~ 600 MeV. If we allow 10% statistical er- 
ror, I think that NLO formula can only apply 
to the data set with Mps ^ 600 MeV, although 
MILC's analysis poses more tight restriction. 

4.2. Wilson-type fermion action 

The CP-PACS collaboration extended their 
coarse (a — 0.22 fm) lattice simulation ^0] to- 
ward smaller quark mass region |18| . The simu- 
lation parameters are tabulated in Tabled The 
lightest quark mass reaches Afy/A/ps ~ 0.35. 
They investigated the chiral behavior of pseudo- 
scalar meson mass, decay constant and PCAC 
quark mass. 

The chiral fit is carried out with the contin- 



uum ChPT formula as well as the Wilson ChPT 
(WChPT) 'SS' formula, which contains the effect 
of scaling violation of 0{a) and 0{a?). The for- 
mer contribution can change the coefficient of the 
chiral logarithm; the latter introduces more sin- 
gular logarithmic behavior and important to re- 
alize the Aoki phase (S^I. Since the NLO calcu- 
lation is not available for the partially quenched 
WChPT with the 0{a^) effect, the ffi is done on 
the unitary data set. Figure[3shows the fits of the 
pseudo-scalar meson mass using the NLO ChPT 
with and without the scaling violation effects. 
The continuum NLO ChPT can fit data only be- 
low Mps ~ 630 MeV (top panel), while the NLO 
WChPT can fit data up to Mpg - 1 GeV. 

The crucial test of the WChPT is to see 
whether it can explain the lattice spacing de- 
pendence of the data, which has been investi- 
gated using the A^^ = 2 data from CP-PACS 
at four lattice spacings From the fits of 

available data at heavier mass region the scaling 
violation is not well described by the WChPT 
formula. They need simultaneous fit on PCAC 
quark masses, decay constants and meson masses, 
and smaller quark mass data including partially 
quenched data at finner lattice spacings in order 
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Figure 5. The continuum ChPT fit to the 
data Mps/Mv < 0.6 (top panel), and the re- 
summed Wilson ChPT fit to the whole data (bot- 
tom panel) ^ . 

to identify the 0{a) and 0{a^) effects in the chiral 
fits. 

The qq-fq collaboration studied the chiral fit 
with the unimproved Wilson action "23I2.41 using 
the NLO WChPT formula including the 0{a) ef- 
fects The simulation parameters are shown 
in Table ^ which cover the similar region as in 
|18| . They also estimate the (continuum) analytic 
NNLO contribution and the low energy constants 
of ChPT. 

Using the double ratio method 4 they esti- 
mate the contribution of NNLO and 0{a) effect. 
Their data indicate that the 0{a) contribution is 
not important, but the NNLO is. ft seems that 
this contradicts with the naive order counting for 
a ~ 0.2 fm (see also 0). Because the qq-l-q data 
point /3 = 5.1 is close to the recently observed 
first order phase transition point f3 = 5.2 |61I62| . 



and one may suspect that the data could be badly 
distorted. Therefore, one needs a detailed study 
of scaling to draw definite conclusion. 

To further discuss this point, we show in Fig- 
ureElthe double ratio against the ratio of valence 
and sea quark masses ^ — mv/ms from the par- 
tially quenched data at a 0.09 fm of the JLQCD 
collaboration [20>63j . RRf and RRn are defined 
as 



RRf 
RRn 



fvs 
fssfvv 
4^ 



(1 



(3) 
(4) 



with V denoting mvai ^ 'm'sca. f^nd S means 
the unitary point 23 . The lines are the the- 
oretically expected chiral logarithm at NLO in 
the continuum ChPT. Filled symbols {MsS ~ 
600 MeV) below ^ = 1 are close to the cor- 
responding theoretical expectation, compared to 
the heavier quark mass data, ft suggests that the 
ChPT is vahd only below Mps < 0.6 GeV. The 
detailed scaling study and application of partially 
quenched WChPT formula including 0{a^ ) effect 
is, however, needed again. 

The UKQCD collaboration added a new data 
at Afps/Mv ^ 0.44 (Mps ~ 400 MeV) [El to 
the previous simulations Large finite vol- 

ume effect is expected, since their physical vol- 
ume is L ~ 1.5 fm and Mp^L — 3.2 at the light- 
est data. They observed some curvature for the 
pseudo-scalar meson decay constant as a function 
of pseudo-scalar meson mass squared as shown 
in Figure \7\ The bursts are the experimental 
values of and /k and the solid curve is the 
continuum ChPT curve fitted to the experimen- 
tal values. The open boxes are the lattice data. 
The crosses are the decay constants in the finite 
volume {L = 1.5 fm) expected from continuum 
ChPT formula They argue that the cur- 

vature observed in lattice data is the result of fi- 
nite volume effect from the chiral logarithm. Fur- 
ther clarification is, however, needed through the 
systematic study on the volume dependence and 
quark mass dependence. 
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Figure 6. RRn and RRf as a function of ^, data from JLQCD, 20^ x 48, Mps =0.6-1.4 GeV HOI 
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12 3 4 

Figure 7. The pseudo-scalar decay constant as a 
function of the pseudo-scalar meson mass squared 
from UKQCD 

4.3. twisted mass QCD and a surprise of a 
first-order phase transition 

Since the chirally twisted mass introduces a 
lower bound on the eigenvalue of the (hcrmitian) 



Wilson-Dirac operator, the unquenched simula- 
tion could be substantially faster for small quark 
masses compared to the (untwisted) Wilson-type 
quarks. Two-flavor twisted mass Wilson quark 
simulations have been started and the flrst results 
are reported by [SJj- Surprisingly, they found a 
strong metastability in the plaquette expectation 
value at /3 = 5.2 for twisted mass of /i = 0-0.1. 
The quark mass determined through the PC AC 
relation changes its sign at the phase gap, and 
the pion mass does not vanish by tuning hop- 
ping parameter. The origin of the metastabil- 
ity may be attributed the possible scenario de- 
scribed in inn], that is that the negative coef- 
ficient of the O(a^) term in the chiral effective 
Lagrangian could cause the absence of the Aoki 
phase. This observation raises a new question 
on the phase structure in the (/?, k) plane in un- 
quenched simulations with Wilson type quarks. 
To reach the chiral limit at finite lattice spacings 
one needs the Aoki phase separated by a second 
order phase transition, but it is not necessarily 
the case. Then, the chiral extrapolation must 
be done after taking the continuum limit. The 
extension of the phase diagram to twisted mass 
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direction is also investigated in the effective La- 
grangian 

4.4. Domain-wall fermion 

The RBC collaboration performed two-flavor 
dynamical simulations with the domain-wall 
quark action '31'22' . It is expected that the 
Ginsparg- Wilson (GW) type fermion has much 
better chiral property because of the exact chi- 
ral symmetry at finite lattice spacings. The do- 
main wall fermion, however, in finite extent in 
the fifth dimension does not obey the GW rela- 
tion exactly. The violation of the symmetry can 
be measured by the residual mass defined through 
the Ward-Takahashi identity of the DW fermion 
and the residual mass remains even in the mass- 
less limit. The lattice size in the fifth dimension 
is A^s = 12, with which the residual mass is small 
enough arrircs — 0.001. Using three quark masses 
rUsea ~ Ws/2, 3ms/4, m^, they analyzed the 
pseudo-scalar meson mass and the decay constant 
using continuum NLO ChPT formula. They ob- 
served that for the pseudo-scalar meson mass the 
partially quenched NLO ChPT can fit the data 
for < 3ms /4 (Mps ^ 630 MeV), while gets 
much worse if one includes the heavier mass data 
(Afps ^ 690 MeV). For the decay constant, the 
NLO formula cannot describe the data well even 
if one restrict the masses in < 3ms /4 (Mps < 630 
MeV) . This is probably because of limited avail- 
able data points compared to that of meson mass 
case. Although they observed f-K/rUp — 0.170(8), 
and //f/m^ = 1.18(1) from the LO fit in fair 
agreement with experiment, the full NLO anal- 
ysis with more data point would be needed to 
control the chiral extrapolation. 

4.5. Summary 

For the KS-type fermions the chiral logarithm 
is observed by the MILC collaboration for small 
sea quark masses using the SChPT formula. On 
the other hand, there is no definite conclusion 
for the Wilson-type fermions primarily because 
the unquenched simulations are still limited to 
relatively heavy quark masses. As a result, chiral 
extrapolation with appropriate ChPT formula in 
finite lattice spacing is not well investigated so 
far. Nevertheless, from the observation from the 



MILC results and the studies discussed in this 
section, it is likely that ChPT can only be applied 
at Afps ^ 0.6 GcV. For Wilson-type fermions, 
development of simulation algorithms that allow 
to enter this region [35] is crucial for progress. 

5. Other topics on hadron spectrum 

5.1. Mixed action simulations 

The Ginsparg- Wilson fermions improve the chi- 
ral property. However, it is numerically so de- 
manding that the unquenched simulations are 
hard task especially for light quarks. Use of differ- 
ent quark actions for valence and sea quarks can 
provide a clue to avoid too large computational 
cost while partly keeping the good chiral prop- 
erty. The partially quenched quark mass combi- 
nations can also be studied in this setup. A pos- 
sible problem is the violation of unitarity, which 
can be removed only in the continuum limit. 

As I already mentioned, the unquenched con- 
figurations with KS-type sea quarks reaches the 
chiral regime. Various measurements on these 
configurations have been carried out. The light 
hadron spectrum in Section [3 . 31 mav also be con- 
sidered as a mixed action calculation because the 
treatment of the KS-Dirac operator is different 
between valence and sea; while sea quarks are ex- 
pressed by the fourth-root trick, valence quarks 
are treated by the quantum number projection. 

Preliminary results for the light hadron spec- 
trum with the overlap or domain wall valence 
quarks on the Nf = 2 + 1 KS-type sea quarks 
generated by the MILC collaboration have been 
reported in ^67.68^ . The valence mass is tuned by 
matching pseudo-scalar meson masses calculated 
with the domain wall and KS quarks |67| . 

5.2. Flavor singlet meson mass 

The large splitting between the flavor singlet 77' 
mass and other flavor octet meson masses should 
be explained by lattice QCD simulations. There 
have been a series of studies calculating the flavor 
singlet mass in unquenched simulations 69 7 01711 

The continuum limit of the 77' mass in Nf = 2 
was studied for the RC(T) action [7J, and M^/ = 
960(87) dais) MeV was obtained where the sys- 
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Figure 8. The rj' meson mass in Nf = 2 as a 
function of Alps /My quoted from Results 
are obtained with PC [TW7n| (crosses), PW [721 
(bursts), PW (squares), and RC(TP) (di- 
amond (/3 = 2.1), octagon (/3 — 1.95), fancy plus 
(/3 — 1.8)) actions. 

tematic errors come from the continuum and chi- 
ral extrapolations. This value seems consistent 
with the experimental value M,y = 956 MeV. 
However, the lattice result represents the value in 
the Nf = 2 world for which the expectation for 
the iso-singlet meson mass is 861 MeV (method 
described in |19I7U| ) or 715 MeV (from Witten- 
Veneziano formula). The above result from the 
RC(T) action is consistent within the error. The 
rather large systematic error in the negative di- 
rection comes from continuum extrapolation. It 
seems difficult to obtain a reliable continuum ex- 
trapolation because of large statistical errors even 
if we add some points with different lattice spac- 
ings around a^^ — 2.5 GeV. One would need a 
better method to improve the statistical signal. 
The chiral extrapolation should be carried out for 
both Mpg and Mps as a linear function of Mpg, 
since reliable fit ansatz is not known. 

Figure |H1 shows the Mps/My dependence of 
Mri' compiled in 19,. The data show a similar 
dependence on Mpg/My except for crosses. The 
data at Hghter two crosses |^ {a = 0.1 fm) seems 
to be constant. Squares |69j are obtained with Z2 



noise method, while bursts |72I74| are obtained 
with truncated eigenvalue approximation (TEA) 
method from the same SESAM configuration [27] 
(symbols are slightly shifted horizontally to avoid 
overlapping). They are consistent with each other 
while the TEA has smaller statistical error. The 
value at physical point is different: ^ 290 

MeV with TEA ,72.74 and - 5201:5^ MeV with 
Z2 [nni- Systematic study of chiral extrapolation 
with more data at smaller quark masses will clar- 
ify the difference of these data. The 10-20% scal- 
ing violantion from ©(oAqcd) error is also ex- 
pected with unimproved gauge and quark actions 
even at ~ 2.3 GeV lattice. 

The effect of rj and r/ mixing can be treated by 
the partially quenched strange analysis within the 
Nf = 2 simulations 72 74,70.73^ . The quadratic 
mass matrix in the quark-flavor basis is obtained 
in |7i^l74l7()| on the a"! ~ 2.3 GeV SESAM con- 
figurations and they obtained — 292(31) MeV 
and AI,y — 686(31) MeV after the diagonarilza- 
tion of the mass matrix. While these individual 
numbers are still below the experimental num- 
bers, the splitting is comparable to the experi- 
mental value. 

5.3. Universality of quenched hadron spec- 
trum 

Sometime ago, it was pointed out that the 
values of Mn/A^v at Alps /My — 0.5 disagree 
in the continuum limit between the PKS action 
and PW (RC(TP)) action 75 . The PKS data 
in aAIy < 1.4 were extrapolated with a func- 
tion Co + ci X a^, and the PW (RC(TP)) data 
in aAIy < 1 were extrapolated with cq + ci x a. 

This year reanalysis is done by adding new 
data with various lattice actions generated since 
2000 j76j. The universality is investigated in the 
quenched approximation by checking the consis- 
tency in the continuum limit among the various 
discretization method using the Bayesian fitting 
approach [77| . All the data are simultaneously fit- 
ted with some constraints. The fit functions are 
polynomials which contains leading and higher 
order terms as a function of aMy (for the KS-type 
action only a^" terms are included) with a com- 
mon continuum limit. The coefficients are con- 
strained to stay within reasonable values by the 
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Bayesian prior. Thus, they obtained a reasonable 
value for the consistent continuum limit. The 
continuum value is la below and 4cr below from 
the previous analysis of the PK and PW actions 
respectively. They also investigated the universal 
continuum limit of Afyri and M^r\ with more 
data (lattice actions) and obtained a reasonable 
values. These tests suggest the difficulty to 
estimate the systematic error of continuum ex- 
trapolation from single lattice action and the im- 
portance of the action improvement with many 
lattice spacings. The similar test for universality 
should eventually be done for unquenchcd simu- 
lations. 

6. Conclusions 

The unquenched QCD simulations with dy- 
namical up, down and strange quarks have made 
significant development over the past year. The 
progress has been more notable for the KS-type 
fermions for which an extensive work has been 
made both in pushing the simulations toward 
small quark masses and in analyzing the results. 
Applying staggered ChPT, it was shown that 
their data are consistent with the expected log- 
arithmic chiral behavior once the pseudo scalar 
meson masses are decreased below Afps < 0.6 
GeV. Worries, however, regarding the field theo- 
retic foundation of the fourth-root trick manda- 
tory in N f =2-1-1 simulations with the KS-type 
fermions still need to be clarified. 

Simulations with the Wilson-type fermions 
yielded an encouraging result that the meson 
spectrum moves progressively close to experiment 
as the number of dynamical quarks increases from 
A^/ = to 2 to 2 + 1 . The largest issue is the long 
chiral extrapolation from heavy quark masses in- 
volved in reaching this result; consistency with 
the expected chiral behavior is not established to 
satisfaction. Hence algorithms which allow simu- 
lations with quark masses as light as used in the 
KS-type fermions are deeply needed. As we en- 
countered for nucleon mass in the quenched case, 
cross checks between the KS-type and Wilson- 
type fermions which will be made possible by such 
algorithms are very important for establishing un- 
quenched predictions. 



Looking further ahead, unquenched simula- 
tions with the Ginsparg- Wilson type quark ac- 
tions with the correct chiral behavior are the most 
challenging issue in future lattice QCD. 
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ence. I also thank S. Hashimoto and A Ukawa 
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